On the scaling approach to electron-electron interactions in a chaotic quantum dot 
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A scaling theory is used to study the low energy physics of electron-electron interactions in a 
double quantum dot. We show that the fact that electrons are delocalized over two quantum dots 
does not affect the instability criterion for the description of electron-electron interactions in terms 
of a "universal interaction Hamiltonian" . 
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The statistical distribution of single-particle energy 
levels and wavefunctions in a chaotic quantum dot or 
disordered metal particle is described by random matrix 
theoryiiS The validity of random matrix theory as a sta- 
tistical description of energy levels and wavefunctions fol- 
lows from the existence of a large parameter, the dimen- 
sionless conductance g of the metal grain or the quan- 
tum doti^ (The dimensionless conductance is the ratio 
of the Thouless energy E'x and the mean level spacing 
A.) The same large parameter g allows for a consistent 
and simple description of electron-electron interactions 
in quantum dots and metal grains, by means of the "uni- 
versal interaction Hamiltonian" , which was proposed by 
Kurland, Aleiner, and Altshulen^ (see also Ref. 0). Ac- 
cording to Ref. to leading order in g, the only relevant 
contributions to the interaction Hamiltonian are the ca- 
pacitive charging energy, the long-range exchange inter- 
action, and the "Cooper-channel" interaction, which is 
responsible for the superconducting instability. 

The justification for the "universal interaction Hamil- 
tonian" follows from the statistics of wavefunctions 0q in 
disordered metal grains or chaotic quantum dots. Wave- 
functions determine the matrix elements of the electron- 
electron interaction, 

Vai3jS = j dridr2(f>a{riy(f>f3{r2)* 

xV{ri,r2)Mv2)(t>s{ri). (1) 

The absence of (long-range) wavefunction correlations in 
chaotic quantum dots causes interaction matrix elements 
to be self-averaging. Most averages are zero, except av- 
erages of "diagonal" interaction matrix elements Vaf3-yS 
where the wavefunction indices coincide pair wise. Re- 
placing interaction matrix elements Vap^s by their en- 
semble average {Vap-ys), only the charging energy, ex- 
change coupling, and Cooper channel interaction remain, 
thus leading to the "universal interaction Hamiltonian" . 
Small non-universal corrections to the interaction Hamil- 
tonian follow from residual wavefunction correlations in 
disordered metal grains or quantum dots, which cause 
small fluctuations of the interaction matrix elements 
Vap-yS around their average. Typically, these fluctuations 
are a factor 1/g smaller than the diagonal matrix ele- 
ments. 

Although the off-diagonal interaction matrix elements 
are a factor 1/g smaller than the diagonal elements, they 



are many, and it is legitimate to ask what their role is. 
This question was addressed by Murthy and coworkers 
using a renormalization-group approach, in a series of 
papers^iiS These authors assumed Fermi Liquid inter- 
actions on time scales shorter than h/E^, and used ran- 
dom matrix theory to describe electron dynamics on time 
scales beyond H/Et. Successively integrating out states 
with highest energy, they found that the "universal inter- 
action Hamiltonian" is stable for repulsive Fermi-liquid 
interactions and for weak attractive Fermi-liquid interac- 
tions, whereas an instability occurs when the attraction is 
sufficiently strong. Remarkably, Murthy et al. found that 
the critical attraction strength is a factor 2 In 2 smaller 
than the attraction strength corresponding to the Pomer- 
anchuk instability in the bulk Fermi liquid, thus creat- 
ing a parameter regime where the bulk system is stable, 
whereas the finite-sized system is nottSi^ 

A renormalization-group treatment of interactions in 
chaotic quantum dots requires knowledge of how (non- 
universal) wavefunction correlations depend on the en- 
ergy difference between the wavefunctions involved. The 
answer to this question depends on the detailed shape 
of the quantum dot and is different for diffusive? and 
ballistioiS electron dynamics (see also Ref. @) . Murthy et 
al. bypass this problem by using the eigenfunction cor- 
relations of a g-dimensional random matrix for all wave- 
functions with energy within Et/2 from the Fermi level, 
treating wavefunctions at larger energies as plane waves. 
Whereas the use of random matrix theory is justified for 
energies far below Et only, it cannot be used to describe 
non-universal wavefunction statistics near the Thouless 
energy. Similarly, residual wavefunction correlations will 
persist for energ ies above Et, which are not accounted 
for in Refs. iGlTl la. A correct treatment of wavefunction 
correlations around Et is important for the renormaliza- 
tion group approach, since most of the renormalization 
of the interaction parameters takes place around that en- 
ergy. 

In this communication we apply the renormalization 
group scheme to the special case of a "double quan- 
tum dot", see Fig. 12 inset. The "double quantum dot" 
consists of two quantum dots of roughly equal size cou- 
pled via a point contact with dimensionless conductance 
g/2^ 1. The Thouless energy Et of the double dot sys- 
tem is equal to gA, where A is the double-dot level spac- 
ing (which is half the single-dot level spacing). The di- 
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mensionless conductances of the two individual quantum 
dots are assumed to be much larger than g, so that ran- 
dom matrix theory and the "universal interaction Hamil- 
tonian" can be used to describe wavefunctions and inter- 
actions in each of the dots separately. The advantage 
of the double dot geometry is that wavefunction correla- 
tions for energy differences near Et can be calculated in 
detail, so that no approximations need to be made upon 
constructing the renormalization group for the electron- 
electron interactions. The analogy between the double 
dot system studied here and the single quantum dot stud- 
ied by Murthy et al. is that in the double dot electrons 
are confined to one dot for times well below H/Et, but 
not for larger tim es, w hereas in the ballistic quantum dot 
studied in R.efs.lfl TTISl thev have a well-defined momentum 
for times below H/Et, but not for longer times. 

Our main finding, to be elaborated below, is that, once 
the correct non-universal wavefunction correlations near 
i?T are taken into account, the instability of the uni- 
versal Hamiltonian occurs at precisely the same interac- 
tion strength as the instability of the double-dot system 
without point contact between the dots. Although this 
conclusion is reached for one specific geometry only, the 
structure of our calculation leads us to expect that the 
same is true for the more general Pomeranchuk-type in- 
stabilities studied by Murthy et al. In other words, we 
expect that the fact that Refs. I6l7l l3 find a parameter 
regime where the bulk Fermi Liquid is stable whereas 
the finite-sized system is not is an artifact of the use of 
random matrix theory to describe wavefunction statistics 
up to a distance Et/'2 from the Fermi level. 

We now describe the details of our calculation. For 
technical convenience, we consider a double quantum 
dot with spinless electrons and with broken time-reversal 
symmetry. Using random matrix theory to describe each 
of the dots separately, the non-interacting part of the 
Hamiltonian for the double quantum dot readsii 

where Hi and H2 are N x N hermitian matrices model- 
ing the Hamiltonians of the quantum dots without point 
contact and H12 is a 2N x 2N hermitian matrix model- 
ing the point contact connecting the two quantum dots. 
The elements of Hi , H2 , and H12 are complex numbers 
taken from independent and identical Gaussian distribu- 
tions. The rows and columns of H are labeled by ro- 
man numbers k = 1, . . . , 2N, where k = 1, . . . , iV and 
k = N + 1, . . . , 2N correspond to the left and right dots, 
respectively. Eigenvalues of H are denoted Sa , the corre- 
sponding eigenvector being written (f>a{k). The size N of 
the random matrices Hi and H2 is of the order of the di- 
mensionless conductance of the individual quantum dots 
and is taken to infinity at the end of the calculation. 
The interaction Hamiltonian has the formic 

Hint = ^Uo{ni + h2f + ^Ui{ni - ^^2)^ (3) 



where hi and 77-2 are operators for the number of elec- 
trons in the two individual quantum dots. The first term 
in Eq. corresponds to a "charging energy" for the 
double dot system, whereas the second term in Eq. © 
represents a dipolar interaction. Upon changing to the 
basis of eigenstates of the non-interacting Hamiltonian 
H, the total Hamiltonian reads 

^ - E ^"^"'/'^ + ^ E y^pysi^l^ii^Js^ (4) 

where "0^ and t})^ are creation and annihilation opera- 
tors for an electron in eigenstate a of the non-interacting 
Hamiltonian H and 

k,l 

where c/j = 1 for k = 1 , . . . , and dfe = — 1 for k — 
N+1,...,2N. 

For large g, wavefunction elements (t>a{k) are indepen- 
dently distributed Gaussian random numbers, 

{(pa{k)(l)p{l)) = ^6kiSf,^, (6) 

corrections to Eq. © being of order 1/g. As a result, in- 
teraction matrix elements are self averaging; fiuctuations 
are of relative order 1/g. Only diagonal elements have a 
nonzero average, 

{Va[3js) — UoSasSpj. (7) 

Replacing the interaction matrix elements by their av- 
erage, we find that interactions are described by the re- 
duced interaction Hamiltonian 

i^int = (8) 

where n — hi -\- fi2 is the total number of electrons in 
the double quantum dot. Notice that, in comparison 
to Eq. ||3J), the dipolar interaction has disappeared be- 
cause the electron wavefunctions are delocalized over of 
the entire double-dot system. Equation ^ is the equiv- 
alent of the "universal interaction Hamiltonian" for the 
double quantum dot; the disappearance of the dipolar 
interaction is the double-dot counterpart of the disap- 
pearance of all non-zero-mode Fermi Liquid interactions 
in the original construction of the "universal interaction 
Hamiltonian" .^-^ 

In order to study the importance of the many resid- 
ual interaction matrix elements that are of order l/g, 
we perform a renormalization group analysis, following 
Refs. 6,7 8 (see also Ref. .13'). This analysis is reminis- 
cent of Anderson's "poor man's" treatment of the Kondo 
problemti^ In order to find the effective interaction for 
electrons at the Fermi level ep: one successively inte- 
grates out states at energies far away from ep. Writ- 
ing the cut-off energy as MA/2, we calculate the change 
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U = 



FIG. 1: Diagrammatic representation of the effective inter- 
action U , to second order in perturbation theory. 



of the effective interaction parameters Uq and Ui upon 
changing M to M' < M within second order perturba- 
tion theory, see Fig. ^ 



Uo{M')-Uo{M) = 



jj2 ^ Y^ npie^,) - npie^) 

k,i tJ.u e^-e"' 
X(/.;(fc)0,,(/)C(O</'.(fc), (9) 



k,l i-iv 

X (/.*(fc)0^(/)0:(O0.(fc). 



(10) 



where np(e) is the Fermi function and the sum over in- 
termediate states is such that only states /i and v with 
at least one of the energies or e,^ in the cut-off region 
M'A/2 < le-ep] < MA/2 are to be included. We omit- 
ted exchange contributions to the effective interaction, 
which are unimportant for the large-g limit we consider 
here. We follow Refs.|^QI3 in replacing the product of the 
eigenfunctions (/)^ and (p^ of the intermediate energies 
and El, in Eqs. (O and ((Tn|l b y its ensemble average. How- 
ever, we deviate from Refs. in keeping the precise 
dependence of the ensemble average of the wavefunctions 



and (f>i, on the energy difference 



Repeating 



the analysis of Ref. |l5j for the Hamiltonian (01, the rel- 
evant wavefunction average is found to be (see also Ref. 

m 



(C(fc)0M(fc)€(O'/'.(O) - 



Sfj.1, + Ski 1 + CTfcCT; 
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(11) 



Substituting this into Eqs. and (|10|l . we find that only 
the non-universal first term on the r.h.s. of Eq. Hll() con- 
tributes in the limit N oo. Changing to dimensionless 
interaction parameters Uj = Uj /A, j = 0, 1, and replac- 
ing the difference equations Q and 11U|I by a differential 
equation, one thus finds 



duQ 
dM 
dui 
dM 



0, 



m ' 



M27r2 



(12) 
(13) 



The flow equations (|12|) and (|13|l are solved with the 



boundary conditions Uj 



C/j/A, j = 0,l, ifM 




FIG. 2: Renormalization group flow of the effective dipolar 
coupling ui as a function of the cutoff M. The solid curve cor- 
responds to the critical dipolar interaction strength ui = —1. 
The other curves are (from bottom to top) for ui = —0.75, 
and 1. Inset: schematic drawing of the double quantum dot. 



oo. Integrating Eqs. H12|) and (|13|l . one finds that uq does 
not fiow, uo{M) ~ uq for all M, whereas 



Ml(M) = 



1 4 Mn 

hi arctan — — 

ui TT 2g 
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4g2 + M 27r2 



2 Mn 

— arctan 

TT g 

-1 



(14) 



In Fig.|51 the solution of the flow equation for the effective 
dipolar interaction ui is shown for various values of the 
unrenormalized interaction ui. 

In order to address the stability of the "universal in- 
teraction Hamiltonian" , we analyze Eq. H14|l in the limit 
M I 0, 



ui{M)= ( ^ + l-^ln2 
ui g 



for M i 0. (15) 



For Ml > —1, the effective interaction strength ui re- 
mains bounded as M | 0. This implies that the corre- 
sponding interaction matrix elements at the Fermi level 
remain of order 1/g, justifying the use of the "universal 
interaction Hamiltonian" for those values of the dipolar 
interaction. It is only for the critical dipolar attraction 
strength ui = —1 that ui diverges upon taking the cut- 
off energy MA to zero. This is precisely at the same 
interaction strength as the location of the instability in 
the absence of inter-dot tunneling. 

The renormalization approach of Mur thy et al. differs 
from ours in two respects. First, in Refs. IdTl l^ there is no 
flow of the interaction parameters for M > g. Second, 
in order to describe the flow for M < g, Murthy et al. 
replace the eigenfunction averag e jni by the average of 
eigenfunctions of a random matrix of size g, 



Sfj.,y + Ski 



(16) 
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m/M^oo) = -l/(21n2) / 

Uj(M^oo) = - 1 



M/g 



FIG. 3: Comparison of exact renormalization group flow 
(dash ed) a nd the flow according the calculation scheme of 
Refs. IdTlSl (solid). Flows are shown for the critical value of 
the dipolar interaction strength mi. 



One then obtains the following flow equations for the 
effective dipolar interaction strength ui: 



dui 
dM 
dui 
dM 



if Af > .g, 

2M^ln2 
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if M < g. 



(17a) 
(17b) 



The flow equations (|17l) agree with the exact flow equa- 
tions for the double dot system only for M <^ g and 
M ^ g, but not for the intermediate range M ^ g. The 
solution of the erroneous flow equations H17|) is 



ui(M) 



Ml if M > g, 

[ur^ + 2(1 - M/g) In 2] if M < g. 

.(18) 

One verifles that in this calculation scheme, the "univer- 
sal interaction Hamiltonian" is stable for ui > —1/2 In 2 
only, so that there is a range of dipolar interaction 
strengths —1 < ui < —1/2 In 2 for which the separate 
dots are stable against the formation of a dipolar charge 



distribution, whereas the coupled dots are not. Our exact 
calculation shows that such a result is incorrect. In Fig. 
13 we compare the flow of Eq. 1)18(1 and the exact flow 
of Eq. I|14|) for the critical value of ui. Although it is 
only in the range M g that the renormalization group 
flow of Refs. loQIS and the exact flow for the double dot 
differ, the flow in the range M g in crucial in determin- 
ing the value of the interactions at which the "universal 
interaction Hamiltonian" becomes unstable. 

Before concluding, we would like to make three re- 
marks about the renormalization group calculation pre- 
sented here. First, the one-loop renormalization-group 
result 117(1 is exact in the large-g limit. This follows from 
the same arguments as used to establish the validity of 
one-loop renormalization group in the work of Murthy 
et al&. Second, in the exact calculation performed here 
all flow of interaction parameters arises from the first, 
non-universal term in the wavefunction correlator 1(11(1 . 
which is off-diagonal in the wavefunction indices. This is 
opposite to the calculation of Refs. 6-7-8, where the flow 
arises from a universal and diagonal wavefunction cor- 
relator. Third, mathematically, the fact that the insta- 
bility of the "universal interaction Hamiltonian" occurs 
precisely at wi = — 1 is a consequence of the Lorentzian 
energy dependence of the flrst, non-universal term in the 
wavefunction correlator l(ll() . A Lorentzian is generic for 
non-universal wavefunction correlations in both diffusive 
and ballistic quantum dotsj^'^" for which and g in Eq. 
((11(1 are replaced by eigenfunctions and eigenvalues of the 
diffusion operator or the Perron-Frobenius operator, re- 
spectively, see, e.g., Ref.^. It is because of this similarity 
that we believe that our calculational scheme, including 
our result for the critical interaction strength, extends to 
the general case. 
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